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ABSTRACT 

This paper deals with a new area of radio oceanography, namely, with 
the determination of the sea surface/geoid and its relative variation. 
Results of the altimeter experiment on Skylab to test the Goddard 
S^ace Flight Center (GSFC) geold are discussed. The spaceborne 
altimeter on Skylab revealed for the first time that the sea surface 
of the world's oceans can be "measured" with an accuracy in the 
meter rai^e. Surface variatioas are discussed as they relate to those 
"computed" from satellite orbital dynamics and ground based gravity 
data obtained from the Defense Mapping Agency — Aerospace Center, 
and other organizations. The recent GSFC geoid has been constructed 
from about 400,000 satellite trackit^; data {range, range rate, angles) 
and about 20, 000 ground gravity observations. One of the last experi- 
ments on Skylab was to measure and/or test this geoid over almost 
one orbit. It was found that the computed water surface deviates 
between 5 to 20 m from the measured one. Further outlined are the 
influence of orbital errors on the sea surface and numerical examples 
are given based upon real tracking data. Theoretical orbital error 
analyses cannot shine any new light on these problems because the 
error values involved are in the meter range and the reliability of 
theoretical error analyses diminishes at these rather small values. 
Orbital he^ht error estimates have been computed for geodetic type 
satellites and are found to be in the order of 0. 2 to 5 meters. Such 
errors are generally larger for the ^ylab due to its relatively low 
orbit, large area, venting and the fact that only conventional Unified 
S-Band data (rai^e and range rate) were available for orbit computa- 
tions. However, for the single revolution used in the analyses, 
orbit error amounted to only a few meters. Furthermore, orbit 
errors tend to be long wavelength type errors (~2000 km and longer) 
thus short wavelength oceanc^raphic features, i. e. , seamounts, 
trenches are clearly visible in the Skylab altimeter data. 
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I. INTRODUCTION 


The geoid is the equipotentlal surface that a uniform and static ocean would 
conform to under the influence of only tlm Eartli's (atmosphere included) gravi- 
tational and rotational forces, oiccluding all other external forces (e.g. luni- 
solar and planetary gravitation, currents, and circulations, wind stresses, 
atmospheric pressure chaises, and variations in density, temperature and 
salinity. ) The real sea surface deviates from the geoid due to variations in 
sea water temperature and density (salinity) together with Coriolis forces 
acting on moving water masses and other forces. These "deviations" have 
been referred to as the "sea surface topography" [ Mhther, 1973 ] and can 
reach values dT up to 2 meters or more. For exiunple, a study by Stommel 
[ 1965 ] indicates the possible existence of a quasi-stationary variation of the 
sea surface of up to 1 1/2 m between the latitudes of ±6St Also, there is a 
rise of about 1 1/2 m in the water surface over 70 to 100 km due to the Gulf 
Stream to quote another example. 

Sea surface variations of this magnitude will be detectable in the future assum- 
ing a further improvement to commensurate accuracies of the geoid, the 
earth's gravity field and our orbit determination capability. It is clear that 
satellite altimetry, if accurate to about 2 to 1C cm can significantly contribute 
to an explanation of the discrepancies which exist, for instance, between 
equipotential surfaces of the Earth's gravity field as obtained from geodetic 
leveling and the determination of the mean sea level using tide gauges. At the 
present time there is ccmtroversy on whether the discrepancies between meas- 
urements of equipotential surfaces obtained from geodetic leveling and mean 
sea level are real or not. Systematic errors (e.g. , ccmsistent sunshine from 
the equatorial regions biasing leveling, or variations from currents) could 
account for the discrepancies [^el, 1975, personal commimication ] . Satel- 
lite altimetry also provides an independent means of decoupling the n<m-tidal 
gravitational forces from the others acting on the ocean surface [Mather, 1973] . 

In addition, a complete new set of precision satellite tracking data will be 
obtained in the form of height measurements by an altimeter spacecraft. This 
means a t^ to thirty fold increase in the number of tracking observaticms, 
which will improve ccmsiderably our knowledge of the Earth's gravity field. 

In principle, there is no difference between ranging data obtained from a 
ground based tracking system and that from an altimeter. Such experiments 
were first conducted during the SKYLAB mission in 1973-1974. 

NASA launched an altimeter equipped spacecraft, namely, GEC^C (r ow GEOS- 
3) [ Vonbun, 1971 ] in April 1975 and more sophisticated *>ceanographic space- 
craft, namely, SEASAT-A is planned for a 1978 launch. These latter two 
spacecraft are part of NASA’s Earth and Ocean Physics Program - EOPAP 
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[ NASA, 1971; Vonbun, 1972 ] . A major part of this effort will be directed to 
the study of the sea surface topography as well as the surface wind fields and 
ocean waves, their direction and height. 

II. THE "MEASURED" SEA SURFACE 

The last Skylab missicm, SL-4 flown from November 1973 until early J'ebruary 
1974, gave a real opportunity to "measure" or partially test the G8FC QEM-6 
Detailed Geoid (ocean surface) over part of one oit>it. Previous to this, geoid 
comparisons were performed only among different models (Goddard Space 
Flight Center, Smithsonian Astn^hysical Observatory and Ohio State Univer- 
sity). Height variations between the different geoidal models as large as 20 m 
were observed together with longitudinal and latitudinal shifts of 500 to 1000 km 
of certain long wavelength features. As will be explained later, orbital error 
problems as well as the overall system accuracy do not permit the detection oi 
ocean surface variations in the one meter region at this time. However, this 
first altimeter data is quite valuable in that it provides the first Independent 
data for comparisons with geoids in ocean areas derived with surface i^nd 
satellite gravity data. 

A. THE GODDARD SPACE FLIGHT CENTER GEOID 

During the past few years, a rather extensive effort was undertaken at 
Goddard to improve the knowledge of the Earth's gravitational field. At 
present there exist six gravitational field models at GSFC. These are 
designated as Goddard J^rth Afcdels GEM 1 through GEM 6. The odd 
numbers designate gravity fields purely derived from spacecraft tracking 
data (Optical, Doppler, Radar, Range and Range Rate, Lasers) and the 
even numbered ones are based upon a combination of satellite tracking 
data and ground based gravity measurements. Of specific n J^reat for 
tills discussion is the recent Goddard Earth Model GEM-C s * \:h, et al. , 

19741 which is based upon a combination of about 400,000 > i'^ilsion track- 
ing data of 27 satellites and 1654 5 ° equal area mean gravity anomalies. 

This model is complete to degree and order 16 with some additional terms 
up to degree 22. Marsh and Vincent [1974 ] have developed a model of a 
rather detailed gravimetric geoid using this GEM-6 gravity field as a 
base (Figure 1). They have computed a detailed gravimetric geoid by com- 
bining the GEM-6 with surface gravity data. About 24, 000 1® x 1® mean 
free air gravity anomaly values c^tained from the Defense Mapping Agency/ 
Aerospace Center were used in the computation. This gravity collection 
was further augmented with data from the National Oceanic and Atmospheric 
Administration (NOAA) and many other sources. The accuracy of this 
geoid (approximately the mean sea surface) is on the order of 3 to 5 meteie 
over most of the northern ocean areas and 20 to 15 meters in areas of the 
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Houthc'rn heinlsphoru where usually no surface gravity data are available 
( Marsh and Vincent, 1U74 I. 

B. THK SKVI.AB "TKST" OF TIIK CIKOID 

In the past, there was n‘ully no economic way to test a global geold via 
metric measurements. Skylab, carrying for the first time an orbiting 
active radar altimeter, made such on experiment possible. Toward the 
end of Skylab fV the astronauts turned on the radar altimeter over almost 
a full orbit os shown In l*'igure 1. Figure 2 depicts both the computed sea 
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surface and the measured one. The computed sea surface refers to the 
departures of the gravimetric geold from a reference ellipsoid (ae * 
C378142 m, 1/’ =• 298.25), The r. m.s. deviation between the computed 
and measured sea surface Is 8 m. These differences reflect the effects 
of three main error sources: 1) orbit errors; 2) errors in the gravimetric 
geold; and 3) errors in the altimeter system. Orbit determination for 
SKYLAB is eeitainly not the easiest one to perform to a high degree of 
accuracy. Relatively large drag and maneuvering (control rocket firing) 
make a precision determination of the orbit impossible as compared to 
other more "geodetic" type spacecraft (Vonbun, 19701. Orbit errors are 
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eetimatod to be in the meter range for this particular case, due to the 
intensive tracking by the Unified S>Band System and the fact that the arc 
U ngih was restricted to a single revolution. Thus, the differences are 
primarily due to the latter two error sources. The good agreement shown 
in Figure 2 is signifloatt proof of the value of a radar altimeter for sea 
surface determination. 

A closer exair.inatlon shows that radar runs 202 and 203 deviate as much as 
10 to 15 m from the computed geoid or sea surface in this case. This 
rather large deviation is primarily attributed to error in the gravimetric 
geoid in this a. ja. Unless one considers deviations smaller than 1. 5 to 
2 m one does not need to make specific reference to the sea surface topog- 
raphy as mentioned earlier. Rirther tidal effects and atmospheric condi- 
tions are not taken into account for the same reason. These factors will 
however play an Important role in the future when overall accuracies will 
reach the one meter level. It should be borne in mind that this is the first 
test of satellite altimetry. 

Even though Figure 2 shows the overall picture, some details seem to be 
significant. For instance, the Marianas Trench is clearly visible on 
altimeter run #212 even on this global scale. 

C. TRENCHES AND SEAMOUNTS 

Looking somewhat closer, one can distinguish such features as ocean 
trenches and seamounts as was shown by McGoogan [1974). Figure 3 
depicts the Puerto Rican Trench and Figure 4, the filarlanas Trench 
(#212) in a very clear fashion. Note that the Goddard geoid does not 
follow these short wave length features. This was never anticipated since 
the smallest wavelength of the GEM-6 model is in the order of 2500 km or 
approximately 22 1/2 degrees on the earth surface. As can be seen, the 
altimeter geoid is about 30 m displaced from the GEM-6 geoid. This is 
due to a combination of orbital height errors, long wavelength errors in 
the GEM-6 geoid and errors in the altimeter system and has nothing to do 
with the trench determination itself. As a matter of fact, this demon- 
strates quite clearly that orbital altimetry is a very powerful tool for 
studying details of the sea surface variatitMis independent of the orbital 
constraints and/or errors. 

This holds at least over distances of, say, 200 to 400 km. It can safely 
be assumed that the spacecraft orbit will not follow such sudden changes 
as these features would require, due to the ’'attenuation" of the gravity 
field errors at spacecraft altitudes. Figure 6 shovns the opposite, namely. 
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a seiimount at the Cupc VoHe iHlrjids. A detailed analysis of the altimeter 
footprint noise characteristics and ground track Indicated that the altim- 
eter was not over land during this particular pans. Therefore, the 
"water hill" shown In the altitude measurement at 1S''41'"30' Is not due to 
a partial sensing of one of the ^ape Verde Islands. These short wave- 
length oceanographic features ■, dlls and valleys) which arc clearly observ- 
able In I-'igurcB 3, 4, and 5 arc above the radar noise and arc a direct 
result of the deflections of the local gravity vector. These deflections are 



Figure 3. Skylab Altimeter Measurements 

caused both by mass variatiomi (shape of the bottom topography, 1. e. , 
subsurface mountain ranges and trenches) and by density variations within 
the Earth. With the radar height noise of approximately one meter one 
certainly can expect to distinguish even much smaller features than those 
shown in these graphs. It is anticipated that many more interesting fea- 
tures of the ocean surface will be discovered from GEOS-3 drta. 

IIL ORBITAL ERROUS AND SEA SURFACE TOIXXIRAPHY 

As mentioned briefly before, if altimetry is to be expected to operate very 
accurately (dm-range or better) on a global basis, the orbit of the altimeter 
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WI 14 4’ 229' 

LONG 141 7* MAP NO 212 149 3* 


ngure 4. Skylab IV Pass 97 Mode 1, 31 Jan. 1974 

spacecraft has to be known to within an equivalent error. Additional conditions 
such as requiring the water masses of the total oceans to be constant over one 
or two days may be necessary so that the sea surface topography can be more 
accurately devc rmined than the orbital height. A gain of a factor of 3 to 4 may 
be possible according to Mather (19751. This could mean that a sea surface 
topography error of say 10 cm could be achieved even though the height of the 
spacecraft may not be known to 30 to 40 cm. If this can be accomplished it 
would certainly be of great help but we should, at the moment, assume a one- 
to-one relationship only to be conservative. Thus, the errors in the orbit 
have to be reduced considerably if the sea surface topography error goal of 2 
to 10 cm is to be obtained in the future. 

Where do we stand at this time as far as orbital uncertainties are concerned? 
Error studies arc certainly important and arc being used in our investigation 
but extreme caution should be exercised in their interpretation since it is very 
difficult to compare "computed" errors with real ones [Vonbun, 19701. For 
the consideration at hand we used actual tracking data and our best orbit com- 
putation model. Figure 6 shows orbital imcertaintics obtained for GEOS-2 
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Mfiure <». Satellite Orbital I'ncertainty EHtimates (OUF) 
or Overlaj) Differences in Meters for GE08-II 

usin^ both do{)pler and optical tr:. eking data with different earth gravity fields. 
Using the ’’overlaii method" (Siry and Stewart, l!)69; Vonbun, 1970 | height 
uncertainties between 0.2 to 4,8 meters have actually iieen oi)tained. These 
values give an indication of the ort)ital height uncertainties one can expect at 
present when a two day rather stable orbital arc is computed using our best 
numerical .rbit computation system and Earth models at GSrc. 

rv. CONCLUSIONS 

In conclusion, it can be stated that: a) it was possible to measure the sea 
surface variations (geoid) along a major part of an orbit; b) the deviations of 
the computed and measured geoid agree quite well for this first test; c) details 
such as trenches and seamounts can easily be identified independent of the 
orbit errors; d) the orbital height errors for geodetic type spacecraft at 
present arc of the order of meters; e) these errors will <;onside ral ly be 
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reduced in the next few years, particularly after GEOS-3 data are fully ana- 
lyzed and thus our present knowledge of the earth's gravity field is improved; 
and, therefore, f) ^ture altimeter spacecraft such as SEASAT-A shumld be 
able to determine variations of the sea surface of the order of dm. 
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